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Abstract. I present results on the expected y-ray and synchrotron fluxes due to neutralino annihi- 
lation in the Large Magellanic Cloud based on 1 1]. 



1. INTRODUCTION 

Decades of observational evidence indicate that galaxies are surrounded by massive dark 
matter halos. Of the many candidate dark matter particles that have been proposed, 
the most popular is the neutralino, X- Products of XX annihilation rapidly decay into 
neutrinos and y-rays, as well as electrons and positrons which emit synchrotron radiation 
in local magnetic fields. At a distance of 50.1 kpc [2], and given a wealth of observations 
indicating the existence of significant amounts of dark matter [e.g.,|2j,|3|], the LMC is an 
obvious choice for neutralino detection. 

The flux of both y-ray and synchrotron emission from neutralino annihilation in the 
LMC depends on the square of the density profile of the LMC dark halo. To obtain the 
density profile we fit the rotation velocity HI data from Kim et al. |4|] and carbon star 
data from Alves and Nelson [5]. We considered a range of density profiles, but here I 
discuss only the Navarro, Frenk and White (NFW)Ja,0] profile, and a shallower profile 
- an isothermal sphere with a core as derived in o by fitting the same LMC rotation 
curve data we used. This profile represents a less ideal scenario. Instead of making the 
assumption of a minimum disk, the stellar disk and gas contributions are included in 
the model, hence the contribution from the dark matter is confined to smaller values 
compared with the minimum disk assumption. More details on the various fits to the 
observed rotation curve are given in (H). 



2. THE y-RAY EMISSION 

Neutralino annihilation proceeds through a number of channels, many of which produce 
y-rays in the final state (see e.g., ||8|]). We considered only the continuum y-ray emission 
due to the decay of neutral pions produced in hadronic jets. Given the uncertainties in the 
dark matter density distribution, we used the approximate Hill spectrum |9|] to calculate 
the number of photons above a certain energy threshold, Ny(E > E t h). The simplified 
spectra derived using this approximation vary from the fit obtained in ifioll by no more 
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FIGURE 1. Left panel: Specific intensity as a function of projected distance from the center of the 
LMC for several background contributions and the neutralino signal above 1 GeV. The signal is shown 
for the NFW profile and the isothermal sphere with a core (IS). The background components are labeled 
as follows: G for galactic, EG for extragalactic, CRs for cosmic rays. Right panel: Same as left panel but 
for an energy threshold of 50 GeV. The background components in this case are labeled as follows: H for 
hadronic and E for electronic. 



than a factor of a few for neutralino masses from 10 GeV up to a few TeV. The y-ray 
emission coefficient j is 

j=N r (E>E th ) { -?4 1 p 2 (r), (1) 
m x 

The quantity (ov)y is the thermally averaged cross section times velocity for annihila- 
tion into y-rays, m% is the neutralino mass, and p(r) is the density profile of the halo. The 
specific intensities and fluxes, discussed in what follows, are obtained by appropriate in- 
tegration of the emission coefficient (along lines of sight and over volume, respectively). 
For more details, see 

EGRET detected a flux of (14.4 ±4.7) x 1(T 8 photons (E > 100 MeV) cm" 2 s" 1 
from the LMC lUTll . Using our fit to the NFW profile and an energy threshold of 100 
MeV, the maximum flux produced by a viable SUSY model is ~ 3.3 x 10~ 9 photons 
cm~ 2 s~ ! , corresponding to m x ~ 50 GeV and (ov) r — 2 x 10~ 26 cm 3 s _1 . While being 
consistent with the flux detected by EGRET, even the maximum predicted flux is almost 
two orders of magnitude too low, suggesting that the primary source is cosmic rays. 
In fact the cosmic ray induced y-rays may be an additional background component to 
consider when trying to detect flux from neutralino annihilation. Following a method 
similar to that presented in [12], we calculated this cosmic ray induced background. 

In Figure [T]l show a comparison of the specific intensity of the different backgrounds, 
as well as that for the NFW profile and the isothermal sphere with a core. The left 
panel assumes an energy threshold of 1 GeV which is appropriate for GLAST, while 
the right panel assumes a 50 GeV threshold appropriate for ACTs. For GLAST the 
relevant backgrounds are the galactic and extragalactic diffuse emission JI^l E3] • For 
ACTs the relevant backgrounds are the hadronic and electronic cosmic ray shower 
contributions IU4U15I1 . though it is worth noting that ACTs can reject hadronic showers 
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FIGURE 2. Angular dependence of the signal-to-noise ratio for a NFW dark matter halo profile. For 
GLAST we use A e ff = 10 4 cm 2 and 1 year of observation. For ACTs we use A e ff = 10 9 cm 2 and 1 month 
of observation. 



with high efficiency. Clearly, the cosmic ray background is not dominant at these energy 
thresholds. 

Observationally, the most relevant quantity is the signal-to-noise ratio and its angular 
and energy dependence. In Figure El I present the signal-to-noise for the flux within 
an angle 6, assuming the NFW profile. Results are shown for three different energy 
thresholds with the neutralino parameters for each chosen to optimize the signal. The 
specifications for each instrument are also indicated. For the noise all the relevant 
backgrounds, including the cosmic ray induced emission, were used. From the figure 
it is clear that focusing on the central regions may help to achieve higher signal-to-noise 
ratios. As expected, for the isothermal sphere with a core (not shown here), the results 
are less optimistic. 

In Fig. |3]l show the parts of the SUSY parameter space that will be detectable. The 
points represent SUSY models produced using the DarkSUSY package [16]. Models 
with m x -(ov)y above a given line are accessible to observation by the corresponding 
instrument, while models in the lower region do not yield a detectable flux. For GLAST, 
the solid line is derived for the optimistic case of a flux sensitivity of ~ 10~ 10 photons 
cm~ 2 s , which is the expected point source flux sensitivity of GLAST for energies 
above ~ 1 GeV lfl7ll . This flux sensitivity corresponds to a 5-<7 detection and 1 year 
of on-target observation. If GLAST achieves its expected sensitivity, then it will be 
able to detect the neutralino signal for a significant portion of the parameter space. 
This is particularly true in light of the recently derived limit m x < 500 GeV iflill 
(vertical arrow in Fig. |3j). If no annihilation signal is detected, the corresponding part 
of the parameter space can be ruled out. For an energy threshold equal to 100 MeV 
we take A e ff = 10 4 cm 2 , 1 year of on-target observation and AQlmc — 1.13 x 10~ 2 sr 
for the solid angle relevant for the collection of noise (since the LMC will be resolved 
by GLAST). With the sensitivity for this energy threshold only a small portion of the 
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FIGURE 3. The minimum detectable (<7v)y versus m x for the NFW profile. SUSY models above each 
curve yield a detectable signal for the instrument and observational parameters assumed. The dashed line 
shows detectability by ACTs with an effective area of 10 9 cm 2 which will be achievable only at high 
energy thresholds (~ 1 TeV). The solid line is derived for GLAST using the expected flux sensitivity 



~ 10 10 cm 2 s ' . The dot-dashed line is derived for GLAST for 1 year of observation, an effective area 
of ~ 10 4 cm 2 , and an energy threshold of 100 MeV. 



parameter space will be detectable (dot-dashed line in Fig. |3j). Clearly, the best chance 
for detection is by using moderate energy thresholds, such as the 1 GeV case, for 
which the backgrounds are relatively low but the source photons are still numerous. 
To identify neutralino annihilation as the origin of the observed flux, the spectrum and 
its characteristic features, such as the cutoff at E = m x , may be useful. In addition, 
monochromatic lines produced by neutralino annihilation (e.g., the 77 line at E = m x ) 
can be useful observational signatures, even though these lines are strongly suppressed 
(see, e.g., d). 

The detectability prospects for existing and upcoming ACTs are less optimistic. The 
commonly assumed specifications are A e ff = 10 8 cm 2 , E t ^ = 50 GeV and 100 hours of 
observation. For these specifications no models will be detectable. A large integration 
time (~ 1 month) and effective area (~ 10 9 cm 2 ) (dotted line in Figure© would improve 
the chances for detection. However, such integration times are fairly long for ACT 
observations and such large effective areas for an energy threshold of ~ 50 GeV are 
beyond the goals of existing and upcoming ACTs. Effective areas of order 10 9 cm 2 
are expected to be achieved by ACTs for energies ~ 1 TeV, but the number of 7-rays 
produced by dark matter annihilation with energies > 1 TeV is expected to be zero for 
models where the upper limit on m x < 500 GeV holds. 
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FIGURE 4. The synchrotron flux in Jy as a function of frequency in the range from 19.7 to 8550 MHz. 
The open triangles correspond to LMC data taken over the period from 1959 to 1991; the dashed line is 
a best fit to the data assuming a decomposition into thermal and non-thermal emission. The dotted line 
corresponds to the non-thermal contribution, and the dot-dashed line to the thermal contribution. The flux 
from neutralino annihilation for m x = 50 GeV, = 5 jlG, and (ffv) ~ 2 x 10~ 26 cm 3 s _1 , is given by 
the thick solid line. The remaining line corresponds to = 18.4 fiG, which is the maximum published 
estimate for the total magnetic field. 



3. THE SYNCHROTRON EMISSION 

Neutralino annihilation not only generates neutral pions, but also a comparable number 
of charged pions. In the presence of magnetic fields, the electrons and positrons pro- 
duced in this way generate synchrotron radiation. Using again the Hill spectrum for the 
quark fragmentation, and taking into account synchrotron losses and inverse Compton 
scattering (ICS) we calculated the synchrotron flux F syn in Jy (1 Jy=10 23 ergs cm 2 
s- 1 Hz" 1 ). 

In Figure|4|l show F syn for m % = 50 GeV and (ctv) = 2.2 x 1(T 26 cm 3 s _1 , for the 
NFW dark matter profile. The thick solid line corresponds to = 5/lG; the lighter 
line represents the results for the maximum total magnetic field estimated in literature 
(18.4 jiG [20]). A compilation of radio observations of the LMC from 1959 to 1991 



in the frequency range (19.7 — 8550) MHz is also shown 112 111 . In the more realistic 
case of = 5 /lG, the neutralino induced synchrotron emission clearly may be part 
of the observed flux, but is lower than the total observed flux, especially at higher 
frequencies. In the low frequency region, the signal increases enough to exceed the 
thermal contribution to the observation. 

Clearly, as was the case for the y-ray flux, a significant contribution in the osberved 
radio-emission must be due to CRs. The best hope for distinguishing the synchrotron 
emission generated by cosmic rays from that produced by neutralino annihilation is to 
use the fact that the density profile of cosmic rays differs from the dark matter density 
profile. While the dark matter halo extends significantly beyond the disk, the density of 



cosmic rays is expected to rapidly decrease at large radii. If low radio frequency observa- 
tions of the LMC are made with high angular resolution, the change from cosmic ray to 
neutralino dominance should be apparent as one moves away from the LMC disk. Cur- 
rently there is no telescope that could carry out observations of the LMC at frequencies 
less than about 20 MHz. Such low frequencies are extremely difficult to observe from 
the ground due to ionospheric absorption and scattering. One promising ground-based 
project is the Low Frequency Array (LOFAR) 1E2I1 which should reach frequencies down 
to ~ 10 MHz and a flux sensitivity of a few mJy in 1 hour. However, the site of the instru- 
ment has not yet been decided, so the LMC may not necessarily be observable. To reach 
even lower frequencies, where ionospheric absorption is very intense, space-based in- 
struments are required. The proposed Astronomical Low Frequency Array (ALFA) II23I1 
should reach down to ~ 0.3 MHz. 
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